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Abstract 

We propose a gauged U(1)h horizontal symmetry around TeV scale that is a subgroup of a 



< 

\q SU(3)h horizontal gauge symmetry broken at O(10 14 ) GeV. The breaking generates right-handed 

Majorana neutrino masses through a SU (3)# sextet scalar. A particular Majorana right-handed 
neutrino mass matrix explicitly determines the remnant U(1)h at low energy which only couples to 

Ph 

b — s and u — r in the gauge eigenstate. The dangerous K — K, D — D mixing and B s — > u + u~ are 

<D 

^ kept to be safe because the relevant couplings are suppressed through high powers of small mixing 

angles in the fermion rotation matrix. Our analysis which applies to the general case shows that 

> 

the Tevatron di-muon anomaly can be explained through the B s and Bd mixing while keeping all 
■^j - the other experimental constraints within 90 % C. L. For the B meson decay, the B s — > / u ± r =F is 

00 

O bound. 



the leading leptonic decay channel which is several orders of magnitude below current experimental 
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I. INTRODUCTION 



Horizontal gauge symmetry was proposed as an extension of the SM gauge symmetries 
to unify all families of quarks and leptons [U [2] . Given the three families of quarks and lep- 
tons, SU(3)h is the most natural choice for the horizontal gauge symmetry Interestingly, if 
one assumes all the SM fermions transform under 3 of SU(3)h, the anomaly free condition 
requires three generations of right-handed neutrinos n^ 1 ' 2,3 [3] while the right-handed neu- 
trinos also play important roles in explaining the origin of neutrino masses. Therefore, the 
SU (3) h horizontal gauge symmetry model provides a natural scheme for the seesaw mecha- 
nism [J] generating small masses for light neutrinos [3J. The Majorana neutrino mass term 
for the right-handed neutrinos explicitly breaks the SU(3)h, thus it is often believed that 
the horizontal gauge symmetry should be broken at a very high-energy scale Mr ~ (9(10 14 ) 
GeV. Then it seems impossible to test the SU{3)h gauge interactions in low-energy exper- 
iments. However, it is not always the case as we will show in detail below. Even if some 
subgroup of the SU (3)# remains unbroken, the right-handed neutrinos can still acquire large 
Majorana masses. 

Since n R transform as 3 under SU(3)h, the Majorana neutrino mass term can arise from 
the vacuum expectation value (vev) of an SU(3)h sextet xe, 

n%(xe}ijn 3 R , (1) 

and Mr = (xe)- The light neutrino mass is given by the seesaw mechanism as m v = 
m^(Mij) _1 m£i. In order to explain the the neutrino oscillation data, suitable (xe) and mo 
are required. For mu and the other SM fermion masses, there must exist octet Higgs under 
SU(3)h in order to accommodate the correct mass hierarchy in quarks and leptons. In 
addition, to minimize flavor changing effects induced by the octet Higgses, we employ a 
scenario with additional Higgses and singlet fermions |5] in which the rri£> and quark mass 
matrices or lepton mass matrix are all independent. By taking a suitable gauge choice of 
horizontal symmetry, we always choose the Mr to be a diagonal matrix, 

(a 0^ 



Mr = (X6> 



B 
yO Cj 



(2) 



The (xe) structure explicitly determines the symmetry breaking. The SU(3)h is completely 



broken in a generic vacuum with A ^ B ^ A). However, with a specific vacuum 

of A = B = C for instance, the vacuum (xe) is invariant under a 5*0(3) symmetry and 
the breaking is SU(3)h — > S0(3)h- Being symmetric second rank tensor under SU(3)h, 
the sextet xe transform as x& ~~ >* U T xeU where U = e leaTa and T a is the generator of the 
horizontal symmetry. A general scheme to obtain the unbroken symmetry is derived from 
the condition that if {T, (xe)} — 0, (xe) is invariant under transformation defined by T. 
To illustrate the feature of our proposal, we take a vacuum as C = —B. This vacuum 
(Xs) = diag(A, B, —B) is invariant under the SU (3) generator A 6 as 



T = A P 



A) (A 
1 
1 



(3) 



Consequently, one can identify the unbroken U(1)h gauge symmetry with generator T 1 and 
it can survive to low energy, for instance O(TeV) which may lead to interesting predictions 
in flavor changing neutral current (FCNC) processes. It was also observed that if there 
exist horizontal gauge interactions, CP violation can be realized with only two generations. 
Explicit examples of CP violation due to U(1)h and SU(2)h was discussed in [TJ. If the 
above U(1)h is broken at the low energy, the horizontal gauge boson exchanges can induce 
additional CP violations [1] at low energies through quark and lepton mixings. 

In the last decades, huge experimental efforts had been made in improving the measure- 
ments on CP violation in the B meson system. Very recently, the D0 Collaboration at 
Tevatron has reported a large charge asymmetry in like-sign di-muon A b st in both B s and 
Bd decays with 6.1 fb _1 . 

N++ — N — 

A ^ Ex P) = N++ _ N - = ~ 9 - 57 ± 2 - 51 (stat.) ± 1.46(syst.) x 10" 3 , (5) 



1 This U(l) T has an unitary equivalent representation. By taking a 45° rotation R between 2 nd and 3 rd 
axes in horizontal space, the vev becomes 



R T (X6)R = 



I A 

B | . (4) 
y B 



The U{1) then becomes T' = R^TR = diag(0, 1, —1). We take, throughout this paper, the basis where the 
Majorana mass matrix is diagonal as Mr = diag(A, B, —B) and the unbroken U(1)h generator is given 
by T in Eq. pi. 
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where N ++ (N ) is the event number for bb — > / u + / u + X( / u X). Such a large di-muon 
charge asymmetry has a 3.2 ex deviation from the SM prediction A*,(SM) = (— 2.3±g;|) x 10~ 4 
[5] and many models have been proposed to account for this anomaly [T7H2T)] . The CDF has 
also measured A b sl = 8.0 ± 9.0 ± 6.8 x 10~ 3 [21], using 1.6 fb 1 of data, which has a positive 
value and large uncertainties. Combining the above two results in quadrature (include the 
systematic uncertainty), we have 



A sl ~ -8.5 ± 2.8 x 10~ 3 . (6) 

At the Tevatron both Bd and B s mesons are produced, hence A b sl is related to the charge 
asymmetries a s f in Bd and B s decays by 2 

A b sl = (0.506 ± 0.043)4 + (0.494 ± 0.043K, . (7) 

New physics (NP) contributions in Bd mixing are strictly constrained (we will show it more 
explicitly in the parameter fit later), so only large NP contributions to the B s mixing (com- 
paring to the other meson mixings) are allowed. For the NP contribution, if the mixing 
in the rotation matrix between mass eigenstate and gauge eigenstate is not huge, then one 
would naturally expect the U(1)h that maximizes b — s mixing as in Eq. (|3]). Indeed, for 
a CKM-like rotation matrix, the gauge boson coupling matrix at the tree level in the mass 
basis goes like 

( A 4 A 3 -A ^ 
A 3 -A 2 1 



G 



y-A 1 X'j 

where A is the Wolfenstein parameter [22J around the order of the Cabibbo angle (A ~ 0.1). 
Clearly, the meson mixings between the first two generation are highly suppressed. The NP 
also couples to leptons. However, their contributions to the B meson decay branching ratio 
to electron and muon are highly suppressed (although A should be replaced by some small 



(8) 



2 If the semileptonic b-hadron decays do not involve CP violating phase, then the charge asymmetry is 
directly related to the mixing-induced CP asymmctris in Bd and B s meson oscillations. af s — — (1.7 ± 
9.1 ± 1.5) x 10 -3 [TT] is measured through the time dependence of B® — > [i + D~X and its CP conjugate 
at DO and ag L = — (4.7± 4.6) x 10 -3 [TU] at the B factory. Nevertheless, we do not use those results here 
due to their large uncertainty. 
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mixing of the lepton rotation matrix) 3 . Therefore, we focus on the phenomenology in the 
B meson mixing and decay. 

The paper is organized as follows: in section [XXJ we propose the specific model in which we 



consider in the paper. In Section III we show phenomenological implications of our model 



on flavor physics which has subsection HI A related to meson mixing and subsection HI B 



related to meson decay. Section IV contains our conclusions. 



II. THE MODEL 



The model starts with a gauged SU(3)h model at extremely high energy. By taking all 
the fermions as 3 under SU(3)h- The particle contents under SU(3)h X SU{2) L x U(l)y is 

1 2 

9l:(3,2, -), ujj:(3,l,g), d R : (3, 1, --) 

t L : (3, 2, -1), e R : (3, 1, -2), n R : (3, 1, 0) (10) 

which is exactly vectorial and the SU(3)h is therefore anomaly free symmetry. It is crucial 
to have right-handed neutrino triplets, n 1 ^ 1,2 ' 3 , for anomaly cancellation [3]. The extension 
to the Pati-Salam unification [23] may be straightforward. 

As we have already discussed the sextet breaking in the introduction, here we focus on 
the Yukawa interactions for the other SM fermions and the Dirac neutrino mass matrix. In 
conventional SU(3)h models, in order to break the SU(3)h as well as the SU(2)l x £7(1), 
one usually introduces one H : (1,2,1), four $ 8 : (8,2,1) to generate all the SM fermion 
mass hierarchies 4 . However, the (8, 2, 1) Higgs will induce large FCNC [21] if the Higgs 
is light. To avoid the too large FCNC problem, another proposal is to introduce (8, 1, 0) 

3 Comparing to the other NP coupling matrix in the gauge eigenstate G = diag (1, 1, a) [26], which has a 
coupling matrix 

/ 1 aA 5 (a-2)A 3 \ 

G ~ aA 5 1 (a - 1)A 2 (9) 

\ (a-2)A 3 (a-l)A 2 a ) 

in the mass eigenstate, the lepton decay branching ratio, especially the one to muon is much more sup- 
pressed in our case. 

4 Another possibility is to consider bulk 5*1/(3) horizontal gauge symmetry broken at one boundary brane 
where all fermion masses are generated at that brane. In this case, one do not need color octet Higgses 
to generate SM fermion hierarchies. See Ref. [57j as an example. 
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Higgs[5]. 

(8, 1,0), H: (1,2,1) (i = «, tZ, e, i/) (11) 

In addition, to generate effective Yukawa couplings, a new set of ST7 (2) l singlet fermions is 
introduced 

U L ■ (3, 1, ^), D L : (3, 1, (3, 1, -2), iV L : (3, 1, 0) 

U R : (3,1,^), (3,1, -§),£*: (3,1,-2),^: (3,1,0) . (12) 

These singlet fermions form invariant Dirac masses and act as messengers to generate the 
necessary Yukawa interactions. We take the up-type quark mass matrix as an example. 
Since the octet Higgs is no longer SU{2) L doublet, qiUR^ is forbidden and the up-quark 
Yukawa interactions only arise as 

Ul^Ur + M U U L U R + q L U R H + X u q L u R H (13) 

where qlUrH is universal. After integrating out the heavy fermion fields U l, Ur, the effective 
up-quark Yukawa coupling reduce to 

^(XJ^ + d^M^qiH. (14) 

The same mechanism also applies to the mass generation of down type quarks, charged 
leptons as well as Dirac neutrinos. By assigning the (3>g ) independently, the mixings 
and masses in different fermion sectors are completely independent for each other and one 
can easily accommodate hierarchies and mixings in SM fermions and the Dirac neutrinos. 
This also enables us to choose the Dirac neutrino mass matrix other than nearly-diagonal 
structure. 

After electroweak symmetry breaking, the effective Yukawa coupling of also arises 

as 

(H)M u ~ 1 ur§^u l . (15) 

Then, the $g exchanges induce FCNC's in general. We have checked that they satisfy the 
strongest constraint from K-K mixing, marginally 5 . However, actual effects depends on the 
mass spectrum of the $g and hence we do not discuss them in this paper. 

5 Since M ~ ($ 8 ) = M z ,/g H ~ 50 TeV, the effective coupling here is (if) /M of 0(1(T 3 ). With additional 
propagator suppression due to 1/M 2 , the amplitude is 0(1O~ 15 ) GeV~ 2 . 
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Another consequence is that both up and down quark mass matrices become Hermitian 

m\ = m u ,m\ = m d (16) 

Thus, the CP violation in strong interactions due to quark mass matrices, 
arg{det(m u ) det(md)} is absent at least at the tree level[S]. In addition, the Hermit mass 
matrices also require the rotations Ul, Ur in the mass diagonalization U\m u \Jn to be equal 
Ul = Ur. In this case, the horizontal gauge boson couples to vector currents of quarks and 
leptons. As a consequence, pseudo-scalar bosons like B s or do not decay to a pair of 
leptons. However, this is only the result of our specific choice of mass generation model for 
quarks and leptons. In the following analyses, we assume more generic rotation matrices 
and Ur and Ul are taken independent for each other to estimate the prediction. 



III. PHENOMENOLOGICAL IMPLICATIONS IN FLAVOR PHYSICS 



The horizontal gauge interaction is real but family dependent. After the mass diago- 
nalization, the other flavor violation entries as well as new CP violation can arise. The 
Lagrangian of gauge interactions is 



-C 



H 



9 H q' ' L Z'^ + L^R 

g H q L (ypTVZ) . . 7 V L ^ + L R 



(17) 



where stands for the rotation for left-handed g-type quarks and T is the generator of 
U(1)h interaction given in Eq. (J3|. 

Flavor changing interactions in the SM can only be measured via electroweak charged 
current interactions. Therefore, for the SM fermion rotation matrixes, only the left-handed 
ones get constrained from the CKM matrix V^V^y = Vqkm an d one cannot determine 
even V£ and respectively. The other rotations are completely unknown. For simplicity 
of the discussion here, we will assume that all magnitudes of the left-handed mixings are 
CKM-like but the complex phases are 0(1) and unconstrained right-handed mixings have 
the similar structure. Therefore, we have the mixing matrix in the mass eigenstates as 



(G 



i\u/d 
L/R 



( A 4 A 3 ^ 
A 3 -A 2 1 
-A 1 A 2 



(18) 



This U(1)h gauge interaction maximizes the mixing in between second and third generations. 
Mixing magnitude in B s , B d and K° is at the order of (1 : A 2 : A 6 ). The D° — D° mixing 
is also at the order of A 6 suppression comparing with B s mixing. This U(1)h is consistent 
with the phenomenological constraints among different meson mixings. If one assume the 
lepton doublet and right-handed singlet rotations are the similar to the quark sector 6 , one 
can also compute the leptonic decay of mesons. For instance, B s — > decay partial 

width has a A 4 suppression. 

A. Meson Mixing 

At the energy scale mj, the effective Hamiltonian responsible for neutral meson mixing 
(and in particular B s — B s mixing) through the tree-level exchange of Z' is 

H = Ci L {m h )Ol L + Cl R {m h )Ol R + C^(m 6 )0^ + G» R {m b )&l R , (19) 

where the AF = 2 operators are given by 

°ll = QiYPLQjqa^PLqj, 0% R = qiYPRQjQil^PRQj, 

°lr = Qil^PLqjqaixPRqj, &[ R = qiP L qjCtiPRqj- , (20) 

and the Wilson coefficients at M z > scale are {C % [ R {M Z ') = 0) 

Cl L {M zl ) = ^{Gff CUM*) = %r{G'tf CUMz.) = ^(G'L 3 ) 2 (21) 

where gn is the horizontal gauge coupling at M z > scale and M Z i is the horizontal gauge 
boson mass. 

In order to calculate the B physics observables, one has to take into account the running 
effects of the four operators above. The relation between these four operators at the M z > 
and m?, scale is presented in Appendix A. After one obtains the Wilson coefficients at m b 

6 Within the minimal SU(5) or 50(10) grand unification theory (GUT), both left-handed and the conjugate 
of right-handed states are embedded into the same GUT multiplet. But in the SU(3)h, left-handed states 
and right-handed states both transform as 3. The GUT multiplet will contain both 3 and 3 and the 
horizontal gauge symmetry model is not consistent with the minimal GUT. Therefore, we will not assume 
any correlation among rotations of quarks and leptons. 
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scale, by using the relevant hadronic matrix elements [39] 



( B l\ LL/RR\ B <l) ~ ^ m B q fB q B LL/RR 

(B q \0%\B q ) » 4 m ^44 



g "V&q 

(BMW « \m Bq f B B^ R (22) 
Here we use + m g ) 2 « 1 and assume 5 L ^, ~ 5^ ~ B RR = B b ^ L = B Bq . Then we 



obtain Mf 2 



Ml 2 = (B q \H\B q ) = - l -f 2 Bq m Bq B Bq 



.(23) 



From the discussion above, the flavor off diagonal coupling between the horizontal gauge 
boson Z' and the first two generation quarks are highly suppressed (For the CKM like 
rotation matrix, it is at least A 3 suppressed), so we will neglect the new physics contributions 
to the K — K 7 and D — D mixing. The Z'-b-s and Z'-b-d couplings, on the other hand, 
is either unsuppressed or A suppressed, hence we expect large new physics contributions to 
modify the magnitudes and the phases of M^ s , where M^ s are off-diagonal mixing matrix 



elements in Eq.23 We can parametrize such effects by 

Mfi s = (Mf 2 /S ) SM A d/s A s = \A d/s \e<*. (24) 

The experimental measured observables are summarized as follows, Am d / S and AT d / s 
measures the mass and decay width difference between the heavy and light mass eigenstates 
of the B d / S mesons 9 . agl is the charge asymmetry in semileptonic B d / S decays. j3 d or /3 a 
measure the time-dependent CP violating phases in the hadronic B decay channel B d — > 



7 The constraint from the CP violation in the K — K mixing is, in fact, marginal for A ~ 0.1, gn — 0.02 
and M z > ^ 1 TeV PQ. 

8 The CP violation phase is defined as <f> = Arg(— M12/T12) so we choose T q to be real here. 

9 The experimental uncertainty in the measurements Ar^ = sign(ReAcp)Ar c ;/r c ( = 0.009 ± 0.037. is too 
big for us to consider it. 
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FIG. 1: The Abs(G' fes ) versus Arg(G' 6s ) region plot for M z > = 1 TeV and g H = 0.02. In the left 
panel, we choose Abs{G' bd ) / Khs(G' bs ) = 0.1 and Arg(G' M ) = Arg(G / {)S ) while in the right panel, 
we choose Abs(G' M )/Abs(G' fes ) = 0.05 and Arg(G' M ) = Arg(G" bs ). The magenta, green, and yellow 
region with dashed, dot-dashed and dotted boundary stands for the allowed parameter space at 
90% C. L. for Am s/d , AT S and S^ K /I3 S . 
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FIG. 2: The similar plots as FIG. [T] but combine both the B d and B s experimental constraints. 
The small overlapped region labeled with "All" is the parameter space that fits for all experimental 
constraints at 90% C.L. 

J /if) Ks or B s J/ip cf). They are shifted by the CPV phases in B d or B s mixing. 

Am d/S = Am^|A d/s | , 
AT, = Arf* cos(^ M + <£) , 

° SL Am™ \A d/s \ ' 
StpK = sin (2(3 d + 4> d ) , 
20?* = 2/3 s - # . (25) 

The theoretical inputs are listed in the Table [!} All the decay constants and bag param- 
eters are used from Ref. [37]. Notice that we use the calculations for Am^ s = 2|Mf 2 | SM , 
Am Bs = 2|Mf 2 | 5M and s in Ref. [36] which uses the more recent decay constants and bag 
parameters with much smaller uncertainties. All of the rest SM inputs are either from [S] 
or [6] 

All experimental measurements which are used to compare with our model outputs are 
listed in Table |h| For the like-sign dimuon charge asymmetry, we use A b SL ~ —(8.5 ± 2.8) x 
10 -3 which combine the D0 measurements with the CDF measurements. For /3f xp and T s 
measured by both CDF and D0 [T2l - [T5] . we use the combined results with each measurements 
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TABLE I: The theoretical input parameters [SJ. 



Pd 




m Bd 




0.38 ±0.01 [9J 


0.018 ±0.001 


5.28 GeV 


5.37 GeV 


fB d V&d 


fB a VBs 




fB a 


(216 ± 15) MeV 


(275 ± 13) MeV 


192.8 ±9.9 MeV 


238.8 ±9.5 MeV 


<t>r 




(Am B J SM 


(Am s J SM 


-0.091^1 


(4.2 ±1.4) x 10~ 3 


0.53 ±0.12 ps" 1 


19.30 ±2.2 ps" 1 


(Ar d ) Ex p 


(Ar s ) Ex p 






(2.67±8;g|) x 10- 3 ps^ 1 


0.098 ± 0.024 ps" 1 







in 2.8 fb _1 [14J where do not add the most recent CDF results [15] since it is difficult for us 
to extrapolated their contributions between 2.8 — 5.2 fb _1 . Nevertheless, we find that our 
conclusion in the paper does not change if we use the combined results in Ref. [20] 10 . 

The experimental constraints on the parameter space of our model are presented in Fig [I] 
and[2j The parameter Abs(G' bs ) and Arg(G' bs ) are quite similar as the parameter h 2 s and a s in 
Ref. [22] except for a overall factor with very small phase related to the (M^ 2 ) SM , therefore 
our results here can be used as the model independent analysis after some re-parametrization. 
For illustration, we choose parameter for Bd couplings Abs(G' bd ) = Abs(G' bs )/10, Arg(G' bd ) = 
Arg(G' bs ) in which there are sizable contribution to A b sl from Bd and Abs(G' bd ) = Abs(G" fes )/20, 
Arg(G' bd ) = Aig(G' bs ) in which the a d sl contribution to A h sl is negligible. In contrast to the 
paper [IB], it is clearly that from the upper rig ht plot in FIG [TJ there is no region allowed 
by all experimental constraints within the la. The best fitted region for the phase 
Aig(G' bs ) C (7r/2, 37r/4) are quite consistent with the one found in Ref. [25]. However, since 
the Ref. [25] essentially marginalize over V 12 in the range — 0.25 ps _1 and use the best fit 
points in which Ar s is about 2.5 times larger than the prediction, the goodness of fit in our 

10 For the /3f xp at 90% C. L., we find /3f xp C (0.13,0.74) from Fig. 7 in Ref. 44 since 1 - CL is almost 
linear to Log(L) where L is the likelihood ratio. 

11 We notice that in [18], they use SM prediction directly from [8] which use the old decay constants and bag 
parameters. Nevertheless, we find that our conclusion still holds because the parameter space mentioned 
in PI] are not allowed by (Ato b J Exp and (Ar s ) Exp . In the other paper which contains the model 
independent fit [T3] , they simply neglect all the large uncertainties from the SM prediction which lead to 
the wrong upper bound of A b sl . 
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result here is reduced significantly comparing to the one in Ref. [25 



TABLE II: The experimental data. 



oExp 


(A) Exp 




(Am Bs ) & P (AT,) Ex p 


0.655 ± 0.024 [TO] 


44+ - 17 


0.507 ±0.005 ps" 


1 17.77 ±0.12 ps" 1 0.154±g;gf 4 ps" 1 


A sl 


B Exp (B s e+e- 


) B E *v(B s -> 


-) S Ex p(5 s -> e ± ^) B Ex p(B° -> /u+r^) 


-(8.5 ±2.8) x 10" 3 ps" 


1 < 5.4 x 10~ 5 


< 4.7 x 10~ 8 


< 6.1 x 10" 6 < 3.8 x 10~ 5 



B. Meson Decays 

The new horizontal gauge boson can also mediate meson hadronic decay or leptonic 
decays at tree level. In this session, we choose to discuss the two leading processes, b — > 
sec and b — > s/j respectively to illustrate how meson decays constrain the horizontal 
gauge interaction. The other transitions are always with additional factors of power of A 
suppressions. 

The effective AF = 1 Hamiltonian responsible for neutral meson decay is 

?4 A / =1 = C3Q3 + C5Q5 + C 3 Q 3 + C,Q 5 (26) 
and the leptonic decay like B® s — > 

'His = + CloQlO + CgQg + C10Q10 (27) 

where i, j = e, /z, r as lepton flavor indices. 

In the case of b — > sec transition, SM contribution at tree level is induced via weak 
charged current with a CKM factor Vb c V* c ~ A 2 . Reading from the effective couplings, 
the horizontal gauge boson mediated b —> sec has a factor of Gb s G cc ~ A 2 . The SM and 
horizontal gauge interaction contributions are at the same order of A and the one can simply 
compare their couplings and gauge boson masses to estimate the ratio. As we discuss in 
previous section, the U{1)h is broken at C(TeV) which results in a suppression due to 
(jjhI '9) 4 ( m w/Mz') ~ 1CT 8 . Consequently, the contribution to b — > sec from new horizontal 
gauge boson is completely negligible. In the case of neutral meson mixing, SM leading 
contribution is from box-diagram while the horizontal gauge boson contribution is at tree 

13 



level. Therefore, even if the new horizontal gauge boson is of order TeV, it is still possible 
to change the AM significantly. For decay process, if there exists SM tree level, the above 
argument then always applies. We won't discuss any constraint from such decays. 

Since the horizontal gauge boson also couples to the leptons, there is again tree level 
contribution to the meson leptonic decay. Within the framework of SM, B s pure leptonic 
decays are realized via the electroweak penguin diagrams with Z/^y* — > £ + £~ and therefore, 
there is no lepton flavor violation at all. The constraints on leptonic decay are mostly on 
leptons directly decaying from B meson. If there exists a r in the final state, r decay will 
complicate the search due to the D ± decays. Therefore, the leading constraints are 

Bt(B s < 4.7 x KT 8 

Br(B s e + e~) < 5.4 x KT 5 

Bi(B s e ± fi T ) < 6.1 x 1(T 6 (28) 

In our model, the horizontal gauge boson has maximal coupling to b, s and /i, r and the 
leading leptonic decay constraint is from B s — > /i ± r =F . But as we mentioned, B s — > / u ± r =F 
does not exist in SM physics at leading order and it is only from the new physics contribution. 
Using(0 | 57^756 I B°) = ifsV^B-- one can compute the decay BR as 

where the Wilson coefficients C 9i i = g^tm^ / 'm 2 z , . To estimate the decay BR, we take 
M z > ^ 10 3 GeV, g H ^ 0.02 and substitute m b = 4.7 GeV, m T = 1.777 GeV, f Ba = 230 MeV, 
rriB = 5.3 GeV, tb = 1.6 ps, we obtain 

Bi(B s -> /j+T*) ~ 1.2 x 1(T 9 . (30) 

12 The only relevant search is from CLEO as of B^ —¥ /i ± r T [30J, 

Bi(B d -»■ /i ± r T ) < 3.8 x KT 5 . (31) 

Due to horizontal gauge boson given, the B^ decay partial width has an additional A 2 
suppression so the bound is prediction well below the experimental bound. One can also 

12 The result is based on assumption from left-handed chiral interaction only which is sufficient to estimate 
the maximal value of leptonic decay BR. 
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estimate the B s — y using the above result. The B s — y has a factor of A 4 

suppression then the prediction is about two orders lower than the current experimental 
bound. 

The other possible rare decay which can be induced by the horizontal gauge boson is the 
FCNC decay in top quark, for instance, t — > c/u + / u ± r =F . However, given the large Mz>, the 
three body decay is highly suppressed. 

Even at the top factory like Large Hadron Collider, it is impossible to observe such rare 
decay event. 

IV. CONCLUSION 

The dimuon asymmetry reported by D0 Collaboration is much larger than the SM predic- 
tion which suggests new sources for CP violation. In this paper, we propose the possibility 
to explain such an anomaly through a tree level exchange of a gauged U(1)h horizontal 
symmetry in B meson mixing. The U(1)h horizontal symmetry is a remnant symmetry of 
SU(3)h broken at Mr ~ (9(10 14 ) GeV through a sextet scalar which gives the neutrino 
mass. Such a U(\)h gauge boson only couples to b — s and /i — r in the gauge eigenstate 
which suppresses all other dangerous meson mixings and B meson decays after the flavor 
rotation. For a general flavor rotation matrix we consider there is a parameter region around 
the phase Arg(G" 6 J C (7r/2,37r/4) which fits the data at 90% C. L. For the B decay, the 
dominate enhanced channel is the B s —y /i ± r =F . Nevertheless, such a enhanced decay channel 
is still one order of magnitude smaller than the current experimental bound. 
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Appendix A: RG Running of the AB = 2 Operators 

In the previous section, the Wilson coefficients are given at the Ms scale while to calculate 
the physics processes involving low energy mesons, one will need to calculate relevant Wilson 
coefficients at the low energy scale. The running contains two steps, the first step is from 
Mz> scale to mt where six flavors contribute to the running of a s , the second step is from m t 
to nib, where only five flavors contribute. We summarize the running effects of the relevant 
AF = 2 operators below [39] . 

The operators belonging to the LL/RR, LR sectors read 

QLL = (^W^V^), 

Qrr = {s a ^P L b a )(s^P L bP), 
Qlr = (s a ^P L b a )(s^PRbP). 

Qlr = (s a P L b a )(sPp R b?). (33) 




[r}(nb)]C R R(tit), 
faai(/i&)] [V22(^b)] \C LR (fi t ) 



(34) 



The variables rji is defined as ratio between strong coupling constant at different scales. 
Given the U(1)h gauge boson is around 10 TeV, we have two different running of a s taking 
into account the threshold correction due to top quark. For the running between the B 
physics scale and top quark, we have 

as\^b) 

For the explicit form of n, we list both the expression at the LO (with subscript (0)) and 
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NLO (with subscript (1)) 

v i0) M = vt /23 



ijMfa) = 1.6273(1 -r^f 3 

Vn(vb) = % /23 , 

J0) f s _ 2 , 3/23 -24/23x 

^21 lA*6 J - gl% _ % J > 

^22 Vb) = % 



-24/23 



^(^J = 0.9250?? 5 24/23 + r/ 3/23 (-2.0994 + 1.1744^) , 

r®M = 1.3875(rf /23 -,- 24/23 ), 

r$(fJib) = (-11.7329 + 0.7829r7 5 )r/ 3/23 - r]~ 24/23 (-5.3048 + 16.2548r/ 5 ) , 

rj$(n b ) = (7.9572 - 8.8822^)r/~ 24/23 + 0.92 50r] 26/23 . (36) 

Similarly, for the running between the top quark mass scale and the horizontal gauge 



boson scale, one can replace all the fib, Ht by fit, Hm z , in Eq. (34), where all the r/s are 

v {0) (^ = vt /21 , 

= 1.3707(1 - Ve )4 /21 , 

nSi^b) = o , 

VniVb) = g^e 721 24/21 ) > 

r®M = % 24/21 , 

V$(Hb) = 0.9219r]- 24/21 +% /21 (-2.2194+ 1.2975r/ 6 ) , 

r®W = 1.38 28(, 24/21 -,- 24/21 ), 

= (-10.1463 + 0.8650r/ 6 )r/ 3/21 + % " 24/21 (-6.4603 + 15.7415r/ 6 ) , 

Vw(l*b) = (9.6904 - 10.6122 % )?7 6 24/21 + 0.9219^ 24/21 , (37) 



and 



* - ^mfi ■ (38) 



17 



[1] T. Maehara and T. Yanagida, Prog. Theor. Phys. 60, 822 (1978). T. Maehara and T. Yanagida, 

Prog. Theor. Phys. 61, 1434 (1979). 
[2] F. Wilczek and A. Zee, Phys. Rev. Lett. 42, 421 (1979). 
[3] T. Yanagida, Phys. Rev. D 20, 2986 (1979). 

[4] T. Yanagida, in Proc. of the Workshop on "The Unified Theory and Baryon Number in 
the Universe", Tsukuba, Japan (1979) p. 95; M. Gell-Mann, P. Ramond and R. Slansky, in 
Supergravity, eds. P. van Nieuwenhuizen et al., (North- Holland, 1979), p. 315. P. Minkowski, 
Phys. Lett. B 67 (1977) 421. 



[5] A. Masiero and T. Yanagida, arXiv:hep-ph/9812225 



[6] C. Amsler et al. [Particle Data Group], Phys. Lett. B667, 1 (2008). 



[7] V. M. Abazov et al. [The D0 Collaboration], arXiv: 1005.2757 [hep-ex]. 



[8] A. Lenz and U. Nierste, JHEP 0706, 072 (2007) [arXiv:hep-ph/0612167| . 
[9] M. Bona et al, |arXiv:0906.0953| [hep-ph]. 



[10] E. Barberio et al. [Heavy Flavor Averaging Group], arXiv:0808.1297 



[11] V. M. Abazov et al. [D0 Collaboration], arXiv:0904.3907 



[12] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. Lett. 101, 241801 (2008) arXiv:0802.2255 



D0 Note 5933-CONF, |http : //www-dO . f nal . gov/Run2Physics/WWW/results/prelim/B/ 

IB5871 

[13] T. Aaltonen et al. [CDF Collaboration], CDF Public Note 9458, August 7, 2008, 
http : //www-cdf . f nal .gov/physics/new/bottom/080724 .blessed- tagged_BsJPsiPhi_ 



update_prelim/public_note .pdf 



[14] CDF and D0 Collaborations, CDF note 9787, D0 Note 5928-CONF, http 



|//www-cdf . fnal .go v/physics/new/bottom/090721 . blessed-bet as_combination2 . 8/ 



D0Note5928_CDFNote9787 . pdf . 

[15] L. Oakes (CDF Collaboration), talk at FPCP 2010, May 25-29, Torino, Italy, |http : //agenda. | 



|inf n . it/ getFile . py/ access?contribld=12&resld=0&materialld=slides&conf Id=2635 

[16] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 



[17] A. Dighe, A. Kundu and S. Nandi, Phys. Rev. D 82, 031502 (2010) [arXiv: 1005.4051 



[hep-ph]]. B. A. Dobrescu, P. J. Fox and A. Martin, Phys. Rev. Lett. 105, 041801 (2010) 



18 



I arXiv: 1005.4238 [hep-ph]]. C. H. Chen and G. Faisel, |a7Xiv:1005.4582 [hep-ph]. C. W. Bauer 



and N. D. Dunn, [arXiv: 1006 1629] [hep-ph] . D. Choudhury and D. K. Ghosh, |arXiv: 1006.2171 



[hep-ph]. J. K. Parry, arXiv: 1006.5331 [hep-ph]. S. F. King, arXiv: 1006.5895 [hep-ph] 



[18] N. G. Deshpande, X. G. He and G. Valencia, arXiv:1006.1682 [hep-ph] 



[19] C. H. Chen, C. Q. Geng and W. Wang, |aTXiv:1006.5216| [hep-ph]. 



[20] Y. Bai and A. E. Nelson, |arXiv:1 007.0596 [hep-ph]. 

[21] CDF Collaboration, Measurement of CP asymmetry in semileptonic B decays, Note 9015, 
Oct. 2007. 

[22] L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 

[23] J. C. Pati and A. Salam, Phys. Rev. D 10, 275 (1974) [Erratum-ibid. D 11, 703 (1975)]. 

[24] S. L. Glashow and S. Weinberg, Phys. Rev. D 15, 1958 (1977). 

[25] Z. Ligeti, M. Papucci, G. Perez and J. Zupan, [arXiv:1006.0432| [hep-ph]. 



[26] V. Barger, L. L. Everett, J. Jiang, P. Langacker, T. Liu and C. E. M. Wagner, JHEP 0912, 



048 (2009) |arXiv:0906.3745 [hep-ph]]. V. Barger, L. L. Everett, J. Jiang, P. Langacker, T. Liu 



and C. E. M. Wagner, JHEP 0912, 048 (2009) [arXiv: 0906. 3745 [hep-ph]]. 
[27] G. Cacciapaglia, C. Csaki, J. Galloway, G. Marandella, J. Terning and A. Weiler, JHEP 0804, 



006 (2008) |a rXiv:0709.1714| [hep-ph]]. 
[28] M. Bona et al. [UTfit Collaboration], JHEP 0803, 049 (2008) [arXiv:0707.0"636l [hep-ph]]. 
[29] E. Barberio et al. [Heavy Flavor Averaging Group], arXiv:0808.1297 [hep-ex]. 



[30] J. E. Duboscq, Nucl. Phys. Proc. Suppl. 144, 265 (2005) |arXiv:hep-ex7 0412022 1 . 
[31] G. Ecker and W. Grimus, Z. Phys. C 30, 293 (1986). 

[32] V. Barger, C. W. Chiang, P. Langacker and H. S. Lee, Phys. Lett. B 580, 186 (2004) 
|arXiv:hep-ph/03100"73] . 

[33] J. Charles et al. [CKMfitter Group], Eur. Phys. J. C 41, 1 (2005) [arXiv:hep-ph/0406184] . 
Updated results used in this paper are from |http: / /ckmfitter .in2p3.fr71 "Results as of FPCP 



2006, Vancouver, Canada. 



[34] E. Barberio et al. [The Heavy Flavor Averaging Group], |arXiv:hep-ex /0603003. 

[35] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996) |arXiv:hep- 



|ph/9512 380|. 



[36] C. W. Bauer and N. D. Dunn, arXiv: 1006. 1629 [hep-ph]. 



[37] J. Laiho, E. Lunghi and R. S. Van de Water, Phys. Rev. D 81, 034503 (2010) |arXiv:0910.2928 



19 



[hep-ph]]. For a web version, see www.latticeaverage.org 
[38] A. J. Buras, |arXiv:hep-ph/9806471| 

[39] A. J. Buras, S. Jager and J. Urban, Nucl. Phys. B 605, 600 (2001) | arXiv:hep-ph/0102316| . 



20 



